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Abstract. We discuss the star formation history of the Galaxy, based on the observations of extremely metal-poor stars (EMP) 
in the Galactic halo, to gain an insight into the evolution and structure formation in the early universe. The initial mass function 
(IMF) of EMP stars is derived from the observed fraction of carbon-enhanced EXP (CEMP) stars among the EMP survivors, 
which are thought to originate from the evolution in the close binary systems with mass transfer. Relying upon the theory 
of the evolution of EMP stars and of their binary evolution, we find that stars of metallicity [Fe/H] < —2.5 were formed at 
typical mass of ~ IOMq. The top heavy IMF thus obtained is applied to study the early chemical evolution of the Galaxy. 
We construct the merging history of our Galaxy semi-analytically and derive the metallicity distribution function (MDF) of 
low-mass EMP stars that survive to date with taking into account the contribution of binary systems. It is shown that the 
resultant MDF can well reproduce the observed distribution of EMP survivors, and, in particular, that they almost all stem 
from a less-mass companion in binary systems. We also investigate how first stars affect the MDF of EMP stars. 
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INTRODUCTION 

Extremely metal-poor (EMP) stars in our Galaxy are the 
relics of the early universe and can be a useful probe into 
the structure formation through merging history and into 
the first stars. Thanks to the large-scale survey of metal- 
poor stars in the Galactic halo such as HK survey fl|] and 
Hamburg/ESO (HES) survey|12|], thousands of stars are 
found with [Fe/H] < —2. Such a large number of EMP 
stars enable us to discuss trends in the chemical com- 
position of these stars with use of the follow-up obser- 
vations of high dispersion spectroscopy. Among them, 
one of the known features of EMP stars is that a large 
fraction (> 20 %) of them show the enhancement of car- 
bon, by far larger than the fraction (a few %) that the 
carbon stars occupy among Population I and II stars. 
Carbon-enhancement can be associated with the mass 
transfer from the AGB companion in close binary sys- 
tems as in the case for metal-rich carbon stars such as 
CH- and Ba-stars, attendant also with the enhancement 
of s-process elements. It is natural that we connect the 
origin of carbon-enhanced EMP (CEMP) stars with the 
survivors of low-mass members in the binary system 
whose companions were once on the AGB and experi- 
enced dredge-up of nuclear products. Indeed, most of 
CEMP stars (~ 70 %) show enhancement of s-process 
elements like barium and lead, although others show nor- 
mal abundance of barium. The former is called CEMP-s, 
while the latter CEMP-no or CEMP-nos. 



It is known that the stars of extremely low- and zero- 
metallicity show peculiar evolution, inherent in low- 
and intermediate-mass models. Because of initial smaller 
abundances of CNO elements, the entropy in the hy- 
drogen burning shell is smaller, and hence, the helium 
flash-convection extends through the tail of hydrogen- 
rich layer, differently from the stars of younger gen- 
erations [3]. Then, engulfed hydrogen is mixed inward 
and burns in the middle of convection to ignite a hy- 
drogen shell-flash, and during the decay phase, the sur- 
face convection deepens to carry out CN elements and 
s-process elements to enrich the surface. It is shown that 
this so-called "Helium-Flash Driven Deep Mixing" (He- 
FDDM) occurs at core helium flash on the red giant 
branch or at helium shell flash on the TP- AGB phase for 
[Fe/H] < -2.5 and M < 3.5M flUH]. The low-mass 
members in the binary systems with the primary stars 
in these metallicity and mass range will show enhance- 
ment of these elements and becomes CEMP-s stars after 
they undergo the mass transfer from the AGB compan- 
ions. Suda et al. [7] have applied the models with the 
He-FDDM to interpret the origin of star HEO 107-5240, 
the most iron-deficient star at that time. 

For stars in the mass range of 3.5 < M / Mq < 6, on 
the other hand, hydrogen mixing, and hence, the He- 
FDDM, do not occur but the third dredge-up operates 
to enrich the surface with carbon as the stars of younger 
generations. In this case, however, the enrichment of s- 
process elements, as suggested from the observations is 



not necessarily accompanied [8]. Accordingly, the sec- 
ondary stars with primary star in this range may show 
up as CEMP-nos stars after the mass transfer; the up- 
per limit of the primary mass range is imposed by the 
hot bottom burning in the envelope, which converts the 
dredged carbon into nitrogen. 

In this work, we first explore the initial mass func- 
tion (IMF) in the early universe based on the binary sce- 
nario for the subclasses of CEMP stars, i.e., CEMP-s and 
CEMP-nos stars. Then we Apply the IMF thus obtained 
to figure out the metallicity distribution function (MDF) 
and compare it the observed distribution of EMP stars 
in the Galactic halo. In deriving the MDF, we study the 
chemical evolution of the Galaxy with take into account 
the merging history of minihalos. In this work, we define 
EMP and CEMP stars as those of [Fe/H] < -2.5 and the 
EMP stars with [C/Fe] > 0.5, respectively; CEMP-s and 
CEMP-nos are the subclasses of CEMP, divided accord- 
ing to whether a star has [s/Fe] > 0.5 or not, where [s/Fe] 
denotes the enhancement of s-process elements such as 
Ba and Pb. We exclude stars with the enhancement of 
r-process elements from the group of CEMP-s. 



IMF OF EMP STARS 

The IMF can be determined from the theory of stellar 
evolution and of binary mass transfer once the fraction 
of the number of CEMP stars among EMP stars, i.e., 
~ 10-20 %, and the ratio of the fraction of CEMP- 
s stars to that of CEMP-nos stars, i.e., ~ 0.33 (9[] are 
given. As stated above, the existence of CEMP stars are 
associated with the He-FDDM at RGB or AGB as well 
as the third dredge-up (TDU). For EMP stars, the occur- 
rence of TDU depends on the initial mass and metallicity 
of a star as discussed in this section. On the other hand, 
the production site of s-process elements for EMP stars 
are assumed to take place in the helium flash convec- 
tive zone of primary stars that experience the He-FDDM 
during thermally pulsating AGB phase. We assume the 
standard 13 C pocket not to work under the EMP circum- 
stance, while this assumption may demand in more de- 
tail discussion. Accordingly, CEMP-s stars are consid- 
ered to be the outcome of He-FDDM and/or TDU in low- 
and intermediate-mass stars, while CEMP-nos to experi- 
ence only the TDU in the models of more massive AGB 
stars. The final fates of EMP stars in the various mass 
range are discussed in this conference proceedings [? ]. 
One should note that massive CEMP-nos stars may also 
undergo the hot bottom burning in the convective en- 
velope. Such stars are recognized as nitrogen-rich EMP 
stars lfTol[Tlll . Since these stars are comparable, or even 
exceed, the CEMP-nos stars in number, the above ratio 
of CEMP-s to CEMP-nos is an upper limit. 



The observed fraction of CEMP, and relative frequen- 
cies of CEMP-s and CEMP-nos stars, can be used as the 
constraint on the IMF of EMP stars by assuming that all 
of these class of stars were born as the members of bi- 
nary system whose primary stars have mass in the range 
of ~ 1 - 7M . We assume the IMF of Millar-Scalo type 
lognormal function in the form; 



1 / (log(m/M md )) 2 



(i) 



where M m( j and Am is the medium mass and dispersion 
of the mass function, respectively. From the theory of 
stellar evolution, we find the upper and lower mass limits, 
M up and Mi ow , to the occurrence of He-FDDM and TDU, 
and also the upper and lower limits, d up and di ow , of 
binary separations that the low-mass member can accrete 
the gas from the AGB companions. Then, with relevant 
values of these upper and lower limits, the fraction of the 
stars of these groups are given as follows: 
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where f(P) and n(q) is the period distribution and 
mass-ratio function of the binary, respectively. Here, we 
choose n(q) = \(q = mi/mi) for simplicity and f(P) 
is taken from Duquennoy 111 211 for stars in the Galac- 
tic bulge. We may set the range of the first integral be- 
tween M low = 0.8 Mq and M up — 3.5M Q in which the 
He-FDDM occurs for the fraction of CEMP-s stars, for 
example. The second integral gives the range of binary 
separation that allows the secondary stars to gain the 
amount of carob-rich matter from AGB companion large 
enough to be observed as CEMP stars. The lower bound 
is determined by the radius of the star when it experi- 
ence He-FDDM and/or TDU on the AGB, while the up- 
per bound is determined by the binary separation that al- 
low the secondary to attain the surface carbon enrich- 
ment [C/Fe] > 0.5 through the accretion from the wind 
of AGB stars and despite of the dilution by the surface 
convection on the secondary. 

The predicted ratios of the fraction of CEMP-s, and 
CEMP-nos relative to the fraction of binaries are shown 
as a function of M m( j of the IMF in Figure Q] In order 
to explain the consistent fraction of these group of EMP 
stars, M mt j ~ 1OM is required, although 4 — 6M Q is 
satisfactory if we consider the fraction of CEMP-s stars 
alone. In conclusion, the IMF of the EMP stars are im- 
plied to be quite different from the present-day IMF that 
peaks at a low mass. This result also suggests that the 
most of the observed CEMP stars are the relics of the sec- 
ondary members in the binary systems because the con- 
tribution of low-mass stars, born as single, is very small 
in such top-heavy IMF. 



survivors at metallicity Z is 
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FIGURE 1. Fractions of CEMP-s and CEMP-nos stars 
among EMP stars as a function of the medium mass M mt j, a 
parameter of the IMF in eq. Q}, with the other parameter fixed 
at Am = 0.33. Solid and broken lines plot the proportions of 
CEMP-s and CEMP-nos stars to EMP stars, predicted from the 
binary scenario, respectively. Dotted line denotes the number 
ratio of CEMP-nos to CEMP-s stars. The observed proportion 
of CEMP-s stars (~ 10 ~ 25%) and the ratio between CENP- 
nos and CEMP-s stars (~ 1/3 ~ 1) are bounded by solid and 
dashed horizontal lines and two arrows indicate the ranges 
compatible with the observations. 



EARLY METAL ENRICHMENT OF OUR 
GALAXY 

In this section, we estimate the MDF of the Galaxy and 
compare with the observed MDF among EMP stars by 
applying the IMF of EMP stars derived above. Since the 
metal enrichment of the Galaxy is caused by massive 
stars that explode as Type II supernovae, we have to 
estimate the ejected mass from supernovae in the early 
phase of the Galaxy. Let us suppose that the Galaxy 
of mass M/, has the N(Z) stars, where Z is the average 
metallicity of the Galaxy at a given time. If AN(Z) of 
new stars are born from the gas clouds and the average 
metallicity increases by AZ, the increment in the mass of 
metals by SNe is given by, 

M h AZ= ( Y(m)% t (m)dmAN(Z). (3) 
Jm sn 

where M$n is the lower mass limit of stars that explode 
as Type II SNe, E, t is the time dependent IMF, normalized 
as Jq '% t (m)dm — 1, and Y(m) is the ejected mass of iron 
as a supernova yield from stars of initial mass m and set 
at 0.07 Mq as a average in this study. For simplicity, we 
assume that the supernova ejecta are instantaneously dis- 
tributed in the whole Galaxy and that the chemical com- 
position of the next generation stars are homogeneous. 
With use of the above equation, the number n(Z) of EMP 
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It is worth note that n(Z) is independent on Z because 
£ t (m) is set to be independent of Z in this study. When 
we change the variable for n from Z to [Fe / H]by simply 
setting Z = Z Q 10P e / H ], 



n([Fe/H])=n(Z) 



dZ 



d[Fe/H] 
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(5) 

This gives the theoretical MDF for EMP stars in the 
whole Galactic halo. 

Using the above MDF, we can compare it with the ob- 
served MDF in the Galactic halo, covered by HES sur- 
vey. In order to obtain the partial MDF in the limited 
sample volume by the survey, we assumed that (1) the 
HES sample includes all EMP giants in the Galactic halo 
within the field of view because the limiting magnitude 
of the survey is large enough to reach farthest giants, 
(2) the HES sample includes the same number of EMP 
dwarfs as derived from the total number of dwarfs esti- 
mated from the number giants by considering the rela- 
tive life times and with the stellar luminosity and the flux 
limit of survey taken into account, and (3) the spatial dis- 
tribution of halo EMP stars is isotropic. 




FIGURE 2. Metallicity distribution function of EMP stars. 
Shaded histogram shows distribution of observed EMP stars 
by HES survey [13]. Two straight lines denote the predictions 
from simple one zone model; solid line for the IMF with M mc i = 
10M© and thin broken line for the same IMF as the Galactic 
spheroid. MDF with the top heavy IMF give a reasonable 
fit to the observations except for the metallicity range below 
[Fe/H] < -4. 

Figure [2] shows the comparison with the observed 
MDF. It reveals that the MDF derived from the top heavy 
IMF of M mc j ~ 10M n is compatible with the observed 
one, while the IMF peaked at low-mass significantly 
overestimates the number of observed EMP survivors. 



This is simply because the low-mass peaked IMF pro- 
duces long-lived, low-mass stars much more than the top 
heavy IMF while producing necessary amount of met- 
als. For high-mass peaked IMF, the number of low-mass 
survivors are smaller by a factor of 1,000 than for the 
low-mass peaked IMF. There are two factors. One is the 
small number of low-mass stars born, which decreases 
the numerator in eq. (4). The other is the large metal- 
enrichment from a single generation of stars, which in- 
creases the denominator in the equation. The fact that the 
observed MDF in the Galactic halo is well reproduced for 
[Fe/H] > —4 implies that the early phase of the chemical 
evolution of our Galaxy gives an additional support to the 
top heavy IMF, derived above from the observed ratio of 
CEMP-s and CEMP-nos. It is to be noted, however, that 
the sharp cutoff at [Fe/H] ~ —4 is hardly to be explained 
within this simple one zone model of chemical evolution 
and this problem is discussed in the next section. 

MERGING HISTORY AND FORMATION 
OF EMP STARS 

The observed cutoff of MDF at [Fe/H] < -4 and two 
known stars at [Fe/H] < —5 suggest that they belong to 
the separate population whose origin is different. In par- 
ticular, the former can be taken as evidence of the hi- 
erarchical structure formation process in the early uni- 
verse. In the current understandings of the hierarchical 
clustering in A CDM universe, the typical mass of the 
primordial gas cloud is ~ 1O 5 M , embedded in the dark 
matter mini-halo of mass of ~ 1O 6 M Il4ll5ll . The cur- 
rently observed MDF should reflect the merging history 
of the gas clouds starting from the primordial cloud to- 
tally devoid of metals. During the formation process of 
the Galaxy, the gas clouds are polluted with metals by 
the stellar yields from massive stars. 

With the current knowledge of cold-dark matter sce- 
nario, we can follow the formation process of the Galaxy 
using the merging tree of mini-halos starting with the 
primordial gas clouds. We may construct the merging 
tree program by using the theory of extended Press- 
Schechter approach lfm[T7ll . In developing a merger tree, 
we adopted a method of Somerville & Kolatt [18]. We set 
the halo mass of M/, = 1O 12 M for the final stage and of 
Mi for the minimum mass of star forming halo by assum- 
ing that the virial temperature attains r v; >(M/) = 10 3 A" 
11411 . We also include the effect of chemical evolution 
with the metal-pollution of gas clouds under the assump- 
tion of one-zone approximation. Thus the gas cloud is 
polluted by a single supernova to be enriched with metals 
as much as [Fe/H] = — 3 .5 . In the subsequent generations 
of clouds, stars are formed with polluted gas. If the pol- 
luted mini-haloes are merged with the primordial clouds, 
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FIGURE 3. Metallicity distribution function of EMP stars. 
Thick solid line: Predicted MDF resulting from merger tree 
with the surface pollution taken into account. Thick dashed 
line: Predicted distribution of pristine metallicity. Thin dotted 
line: Distribution of observed EMP stars by HES survey 11311 . 
Predicted MDF is consistent with observation of EMP stars, but 
inconsistent for HMP stars of [Fe/H] < —4. 



the metallicity in the massive mini-halo can be less than - 
3.5 in [Fe / H]. In the primordial clouds of small masses, 
stars are also affected by the accretion from the inter- 
stellar gas polluted with metals due to their low virial 
velocity. We treat the accretion process according to the 
Bondi's formalism with a given virial velocity and the 
stellar mass of O.8M . Depending on the merging his- 
tory of individual clouds, the stars born in the primordial 
mini-haloes can be polluted up to [Fe/H] ~ —5 through 
the accretion of interstellar gas. More details on the pro- 
gram are described at Komiya et al. (2007, in prepara- 
tion). 

The results of the computations are illustrated in Fig- 
ures [3] and [4] Fig. [3] show the results obtained by as- 
suming the IMF of M m( j = 1OM . Clearly, we can see 
the two distinct groups of stars: the lower metallicity 
group comes from the first generation stars, while the 
more metal-rich group comes from the second and subse- 
quent generations. The predicted MDF well reproduces 
the sharp drop at [Fe /H] < —4, although the cutoff may 
is not sharp enough compared with the observation. Seri- 
ous discrepancy is found for the population of primordial 
stars with metal pollution by accretion. Since the merg- 
ing trees starts from a large number of primordial clouds, 
low-mass primordial stars are significantly overproduced 
and survive to date in the Galactic halo. 

Fig. |4] shows the results for different IMF of M mc ( = 
100 Mq. This is the consideration of the suggestion from 
the star formation theory that primordial clouds prefer 
to be born as very massive stars as much as ~ 1OOM . 
This result denoted by thick solid line fits better to 
the observed MDF, althogh the width of the distribu- 
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FIGURE 4. Same as fig. [3] except for stars formed in zero 
metal gas. Thick lines: Case of M m( j = 100 Mm in IMF for zero 
metal gas, and later M m( j = IOMq. Solid and dashed line shows 
MDF after and before surface pollution, respectively. The pre- 
dicted number of Pop. Ill stars is comparable to observed num- 
ber of HMP stars. Thin solid line: All Pop. Ill stars assumed 
to explode as PISNe and all gas to escape from mini-halo. 
The predicted number of EMP stars is much fewer than the 
observed number because PISN eject larger amount of metals. 
Thin dotted line shows the observed MDF. 



tion function for Pop. Ill is still inconsistent. For very 
top heavy IMF, the contribution of pair-instability super- 
nova (PISN) in the mass range of 150 — 28OM becomes 
large. The thin solid line in Fig.|4]is the result of the simu- 
lation if all Pop. Ill are explode as PISNe and eject 10M Q 
of iron per event. Since the velocity of PISN is very large 
enough to exceed the escape velocity of mother clouds, 
the ejected mattar are dispersed into the halo and enrich 
the metallicity of other clouds. As a result, the Pop. Ill 
survivors become metal-rich and merged into the EMP 
population. The slope of the MDF also becomes less 
steep due to the rapid increase of metallicity in the fi- 
nal halo. In conclusion, the massive IMF is prefereble in 
reproducing the observed MDF, but the contribution of 
PISN should be negligible or PISN do not work in the 
early universe. 

SUMMARY 

We discuss about the star formation history of our Galaxy 
using the observation of EMP stars with the current 
knowledge of the stellar evolution, and within the frame- 
work of hierarchical clustering in A CDM universe. 
From the comparison of observed CEMP, CEMP-s, and 
CEMP-nos stars with the binary scenario, the currently 
observed EMP survivors were born in the binary system 
with the IMF of the typical mass of 1OM . We have de- 
rived the metallicity distribution function based on the 
simple one zone model of chemical evolution to find 



that the derived top heavy initial mass function is consis- 
tent with the observed metallicity distribution function. 
We also construct the merging tree of mini-halos starting 
with the primordial gas clouds of ~ 10 5 M Q until the for- 
mation of the present Galactic halo of 1O 12 M . The pre- 
dicted metallicity distribution function with merging his- 
tory can reproduce a sharp cutoff around [Fe/H] ~ —4, 
but is subject to significant overproduction of Pop. Ill ob- 
jects, not supported by the observations. The inclusion of 
the contribution of pair-instability supernovae is not de- 
sirable in reproducing the metallicity distribution func- 
tion of EMP halo stars due to their large amount of iron 
ejecta. 

Although this work is an important step, more sophis- 
tication of the models is desirable for better understand- 
ing of the star formation history of the Galaxy. The next 
generation of the project of the search for EMP stars 
in the Galaxy (SDSS/SEGUE, LAMOST) will provide 
more insight into the understandings of the structure for- 
mation of the Galaxy, leading us to the field of "near-field 
cosmology". 
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